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Abstract 



We study the solvability of the time-independent matrix Schrodinger differential 
equations of the quantum Rabi model and its 2-photon and two-mode generaliza- 
tions in Bargmann Hilbert spaces of entire functions. We show that the Rabi model 
and its 2-photon and two-mode analogs are quasi-exactly solvable. We derive the 



exact, closed- form expressions for the energies and the allowed model parameters 
for all the three cases in the solvable subspaces. Up to a normalization factor, the 
eigenfunctions for these models are given by polynomials whose roots are determined 
by systems of algebraic equations. 
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1 Introduction 



The quantum Rabi model describes the interaction of a two-level atom with a single 
harmonic mode of electromagnetic field. It is perhaps the simplest system for modeling the 
ubiquitous matter-light interactions in modern physics, and has applications in a variety of 
physical fields, including quantum optics [1], cavity and circuit quantum electrodynamics 
[2, 3], solid state semiconductor systems [4] and trapped ions [5]. 

Recently Braak [6] presented a transcendental function defined as an infinite power 
series with coefficients satisfying a three-term recursive relation, and argued that the 
spectrum of the Rabi model is given by the zeros of the transcendental function. This 
theoretical progress has renewed the interest in the Rabi and related models [7, 8, 9, 
10, 11, 12]. However, since Braak's transcendental function is given as an infinite power 
series, unless the model parameters satisfy certain constraints for which the infinite series 
truncates, its exact zeros and therefore closed-form expressions for the energies of the 
Rabi model can not be obtained even for those corresponding to the low-lying spectrum. 
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This is not surprising because, as pointed out in [8, 12], the Rabi model is not exactly 
solvable but quasi-exactly solvable [13, 14, 15]. A typical feature of a quasi-exactly solvable 
system is that only a finite part of the spectrum can be obtained in closed form and the 
remaining part of the spectrum is not algebraically accessible (i.e. can only be determined 
by numerical means). 

In this paper we examine the solvability of the quantum Rabi model and its multi-mode 
and multi-quantum generalizations within the framework of Bargmann Hilbert spaces of 
entire functions [16]. We show that the Rabi model and its 2-photon and two-mode 
generalizations are quasi-exactly solvable. We derive the explicit, closed-form expressions 
for the energies, the allowed model parameters as well as the wavefunctions for all the 
three cases once and for all in terms of the roots of the systems of algebraic equations. We 
note that the energies for the Rabi and 2-photon Rabi models were obtained previously 
[17, 18, 19] by different methods. 

The work is organized as follows. In section 2 we recall the widely accepted charac- 
terization of solvability of a linear differential operator in terms of invariants subspaces. 
In sections 3, 4 and 5, we obtain the exact solutions of the Rabi model and its 2-photon 
and 2-mode generalizations in their respective solvable subspaces. 

2 Exact solvability versus quasi-exact solvability 

Exact solvability and quasi-exact solvability are closely related [20]. A quantum me- 
chanical system is exactly solvable in the Schrodinger picture if all the eigenvalues and 
the corresponding eigenfunctions of the system can be determined exactly, In contrast, 
a system is quasi-exactly solvable if only a finite number of exact eigenvalues and eigen- 
functions can be obtained. Among various characterizations of solvability, the one about 
the existence of invariant polynomial subspaces is conceptually the simplest. 

A linear differential operator H is called quasi-exactly solvable if it has a finite- 
dimensional invariant subspace VV with explicitly-described basis, that is, 

HVm C VV, dim V^ < oo, VV = span{^, • • • , £ d imv^}- 

An immediate consequence of this characterization of quasi-exact solvability for the op- 
erator H is that it can be diagonalized algebraically and exact, closed-form expressions of 
the corresponding spectra can be obtained in the (solvable) subspace VV- The remaining 
part of the spectrum is not analytically accessible and can only be computed through 
approximations (though sometimes rather accurately). If the space VV is a subspace of a 
Bargmann-Hilbert space of entire functions in which H is naturally defined, the solvable 
spectra and the corresponding vectors in VV give the exact eigenvalues and eigenfunctions 
of H, respectively. 



A linear differential operator H is exactly solvable if it preserves an infinite flag of 
finite-dimensional functional spaces, 

Vi C V 2 C • • • C W C • • • , 

whose bases admit explicit analytic forms, that is there exists a sequence of finite- 
dimensional invariant subspaces VV, Af = 1, 2, 3, • • •, 

HV/s C 1/V, dim Vat < oo, W = span{£i, • • • , CdimVjv}- 

As will be seen in next sections, based on the above characterization of solvability of a 
quantum system, the Rabi model and its 2-photon and two-mode analogs are not exactly 
solvable, contrary to the claims by [6, 7, 10, 11]. Rather they are quasi-exactly solvable, 
as was also noted by Moroz [12] for the Rabi model, because only a finite part of their 
spectra can be determined exactly. 

3 Quasi-exact solvability of the quantum Rabi model 

Quasi-exact solvability of the Rabi model has recently been noted by Moroz [12]. Special 
exact spectrum of the model was obtained in [17, 18, 19] by different methods. In this 
section we re-examine this model within the framework of the Bargmann Hilbert space 
of entire functions We will solve the time-independent Schrodinger matrix differential 
equations by means of the functional Bethe ansatz method [21, 22, 23]. In addition to 
the exact spectrum, we are also able to obtain the closed form expressions for the allowed 
model parameters and the polynomial wavefunctions in terms of the roots of a set of 
algebraic equations. 

The Hamiltonian of the Rabi model is 

H = ooa)a + Aa z + ga x o' + a , (3.1) 

where g is the interaction strength, a z , a x are the Pauli matrices describing the two atomic 
levels separated by energy difference 2A, and o' (a) are creation (annihilation) operators 
of a boson mode with frequency ui. In the Bargmann realization a 1 — > z, a — >■ 4-, the 
Hamiltonian becomes a matrix differential operator 

H = ooz- — h Aa z + g a x [ z + z— ) . (3.2) 

az \ az J 

Working in a representation defined by o x diagonal and in terms of the two-component 

(ib (z) \ 
, the time-independent Schrodinger equation gives rise 
i)_(z) j 

to a coupled system of two lst-order differential equations, 

(uz + 9)j-^+{z) + (gz - E)i> + {z) + Ai/>-(z) = 0, 

(cuz - 9)^-(z) - (gz + E)^-{z) + Aij + {z) = 0. (3.3) 



If A = these two equations decouple and reduce to the differential equations of two 
uncoupled displaced harmonic oscillators which can be exactly solved separately [24]. For 
this reason in the following we will concentrate on the A^O case. 
With the substitution ip±(z) = e~ gz ^(f)±(z), it follows, 



fr» + »>sr(£ + B l 



:■"-»>;§- 292 -!j + b 



4> + (z) = -A<j>_(z), 

4>-{z) = -AMz) 



The differential operator 



(3.4) 



C l R = (ujz + g) 



dz 



9_ 

UJ 



E 



(3.5) 



in the first equation is then exactly solvable. Eliminating (f)-(z) from the system we obtain 
the uncoupled differential equation for (f>+(z), 



(ujz-g)-^-- [2gz-^- + E 

dz V u , 



(uJZ + g)±-( 9 - + E 

dz \ uj i 



fa 



(3.6) 



This is a second-order differential equation of Fuchs' type. Explicitly 

d 2 



(ujz- g)(ujz + g)—— + 



dz 2 



-Ibjgz 



w 



2 2g 2 -2Eu)z + -(2g 2 -co 2 ) 



+ 



9' 



2gl^ + E\z + E' 



A 2 - 9 - 

cu 2 



UJ 

(j) + = £0+ 



dz 



0. 



(3.7) 



It is easy to see that for any positive integer n, 



Cz n 



—2ncog + 2g 



K)l 



v n+l 



+ 



n(n - 1)uj 2 + (w 2 - 2g 2 - 2uE)n + E 2 - A 2 - - 



CJ 



+lower order terms. 



(3.8) 



Due to the z n+1 term on the right hand side of (3.8), the operator C is not exactly 
solvable. Only if n — j\f such that —2j\fujg + 2g (E + ^-j = 0, the z n+1 term disap- 
pears from the right hand side, and C preserves a finite dimensional subspace Vat = 
span{l, z, z 2 , ■ ■ ■ , z }. Therefore C is quasi-exactly solvable with invariant subspace 
Vjsf = span{l, z, z 2 , ■ ■ ■ , z^}. 

We now seek exact solutions of the differential equation in the solvable subspace VV- 
Obviously they are polynomials in z of degree jV, which can be written as 



M*) 



IE 



Z - Zi) 



AT =1,2, 



(3.9) 



where Zi are the roots of the polynomial to be determined. The differential equation (3.7) 
is in fact a special case of the general equation studied in [23]. Applying the results of 
this reference, we obtain the energies 

E = u[M- 9 ^ (3.10) 

and the constraint for the system parameters 

A 2 + 2A/V + 2ugJ2zi = Q - ( 3 - n ) 

i=i 

Here Zi satisfy the set of algebraic equations 

* _j_ = 2^? ± (w - iw- r ± g(g - V)/-> j = 1i2i ... iM (3 . 12) 

J^i Zi - Zj [uzi - g) (uzi + g) 

The corresponding wavefunction component i[)+(z) of the model is then given by 

JV 
Mz) = e-" z T[(z-z l ), (3.13) 



i=i 



and the component ip-(z) = e~ 9Z ^ UJ (p^(z) with (f)-(z) determined by the first equation of 
(3.4) for A^O, i.e. (f>-(z) = — ^£jj0 + (z). Because of the exact solvability of C X R) <t>-{z) 
automatically belongs to the same invariant subspace as <f)+(z). 

As an example to the above general expressions, let us consider the M = 1 solution. 
The energy is E — u (l — ^ )• (3.12) becomes 



LLOgZ\ ■+- OJ'Zi -t- 

which has two solutions 



2uugz 2 + u?z x + -(w 2 - 2g 2 ) = 0, (3.14) 



9 2g 2 -u 2 
z\ = , — • (3.15) 

u lug 

Substituting into (3.11) gives the constraints A = and A 2 + 4g 2 = u 2 , respectively. 
The constraint A = corresponds to the case of degenerate atomic levels. The other 
constraint agrees with that obtained in [19] by a different approach. The corresponding 
wave function is 

*♦«=«-*'(*- ^tr) ■ (3 ' 16) 

4 Quasi-exact solvability of the 2-photon quantum 
Rabi model 



The Hamiltonian of the 2-photon Rabi model reads 



,t^2 , 2 



H = Loa^a + Aa z + g a x (a T ) + a 2 . (4.1 



Introduce the operators K±, K by 



K + = \^f, K_ = \a\ K,= l -[a^a+ 1 ^. (4.2) 

Then the Hamiltonian (4.1) becomes 

H = 2u (k - i) + Aa z + 2ga x (K + + #_). (4.3) 

The operators K±, Kq form the usual su(l, 1) Lie algebra. The quadratic Casimir of the 
algebra, C = K + K_ — Kq{Kq — 1), takes the particular values C = ^ in the representation 
(4.2). This is the well-known infinite-dimensional unitary irreducible representation of 
su(l, 1) known as the positive discrete series T> + (q) , where q — \,\- Thus the Fock- 
Hilbert space decomposes into the direct sum of two subspaces 1-L q labeled by q = 1/4, 3/4. 
The basis state in the subspace "H 9 , denoted as \q, n), n — 0, 1, 2, • • •, has the form 

\q,n)= . J ; |0), (4.4) 



2(n + q-\ ] 

and the action of {K±, Kq} in this representation is given by 
K \q,n) = (n + q)\q,n), 
K+\q,n) = J(n + q+-) (n + q + -)\q,n + l), 



K-\q,n) = \l(n + q--J\n + q--J\q,n-l). (4.5) 



Using the Fock-Bargmann correspondence, 



a) — > z, a — > — , In) — > — =, (4.6) 

dz 



we can make the following association 



\q, n) — ► V q , n (z) = (4.7) 

^(2(n + g-i; 



Thus in each subspace % q of the Bargmann space with basis vectors \J/ ?;n (z), the operators 
K±,Kq (4.2) are realized by single- variable 2nd differential operators as 

d z d 2 d , 

K B = z - +q , K + = -, K^2 z - 2+iq -. (4.8) 

In terms of the differential realization, the 2-photon Rabi Hamiltonian becomes 

( d 1\ A (z d 2 d\ tA N 

H = ^ (*& + " ~ 4 j + **' + 29 "' U + ^ + iq ^) (49) 
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Working in a representation defined by a x diagonal and in terms of the two component 
wavefunction, the time-independent Schrodinger equation leads to two coupled 2nd-order 
differential equations, 



d 2 d 

Agz—^z) + (2uz + 8gq)—^ + {z) + 



(jz + 2u;[q--) -E 



i> + (z) + Ai>_(z) 



0. 



d 2 d 

Agz—i>_(z) + (-2uz + 8gq)—^_(z) + 
dz 2 dz 



gz-2u(q 



E 



il>-(z) -Aip + (z) = 0. 
(4.10) 



If A = these two equations decouple and reduce to the differential equations of two un- 
coupled single-mode squeezed harmonic oscillators which can be exactly solved separately 
[24]. For this reason in the following we will concentrate on the A^O case. 
With the substitution 



where 



2£ 



< 1, it follows, 



—?-(l-(l)z i \ 

e 4 (> y ' <p±{z) 



n 



Ag 2 



OJ* 



(4.H) 



z 
d 2 



^9Z— + [2uVtz + 8gq] — + 2quVt - -u- E}ip. 



d o o 1 
.- — h 2quj\l c 

dz ' 2 

d_ w_ 2 
1 dz g 



-A(p. 



Agz— + [2u(n - 2)z + 8gq} — + —(1 - Q)z + 2qu(Vt - 2) + -u + E } <p_ = Ay?+. 
dz 2 dz g 2 J 

(4.12) 



Then the differential operator 



d 2 d \ 

£-2- P = ^gz— + \ 2ooVLz + 8 9Q}j- + 2gwfi --u-E 



(4.13) 



appearing in the first equation is exactly solvable. Eliminating tfi-(z) from the system, 
we obtain the uncoupled differential equation for (p+(z) 

{ d 2 d J 2 1 1 

\ Agz—- + \2u(Q - 2)z + 8gq} — + —(1 - Q)z + 2qu(Q -2) + -u + E} 
[ dz 2 dz g 2 J 

id 2 d 1 1 

x {Agz-— + \2uQz + 8gq]— + 2quQ - -u - E}ip + (z) = -A 2 <p + (z).(4.U) 
dz 2 dz 2 



This is a 4th-order differential equation of Fuchs' type. Explicitly, 



le^dV+eV 



^-l)z 2 + {q+\)z 



dz A 

4u 2 (Q 2 -3Q + l)z 2 + 16cug 



dz 3 



3(q + -)Q-3q-l 



z + 64g 2 q(q+-)} 



d 2 ¥+ 
dz 2 



u 



+ {2—n(i-n)z 2 + 
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8u 2 q{l-Q) + 8u 2 {q + -){l-Q) 2 



+4u{E - 2u{q + -)) 



z + 32ujgq 



(q + -)Q-q 



J dz 



— (1 - ft) [2gwfi - -w - E ) : 



1 



+4w 2 g 2 (1 -nf-(E- 2u(q - ±) ) + A 2 ^+ = 0. 



(4.15) 



By the arguments similar to those in the previous section, this equation is quasi-exactly 
solvable provided that the system parameters A, a; and g satisfy certain constraints, and 
exact solutions are polynomials in z in the solvable sector. We thus seek polynomial 
solutions of the form to the above differential equation, 



<P+(z) 



M 
IK 



Z - Zi 



M = 1,2, 



(4.16) 



where Ai is the degree of the polynomial solution and Zi are the roots of the polynomial 
to be determined. Substituting into (4.15) and dividing on both sides by (f>+(z) yield 

E-2u(q--)\ -A 2 -Au 2 q 2 {l-Q) 2 

M 1 M A 

= ^ 2 z 2 Y.-^- E 



i=l Z Z i V +l+j+i 



CO ,„ „, 2 , 1 



_(n-iy + ( g + -> 



+ 64/ 

+ J4w 2 (fi 2 - 3fi + l)* 2 + I6cug 



. (Zj — 2 p )(2j — 2j)(.Zj — ^j) 
1 



M 

E 



E 



i=i z z i i&# 



(Zi - Z t )(Zi - Zj) 



3(q+-)Q-3q-l 



+64/g(g + i)}E- 1 -E 



i=l ^ Z * j^i Z i Z j 



U 



+ {2—n(i-n)z* 

I 9 



8w 2 g(l - ft) + 8u 2 (q + -)(1 - ft) 2 



+4w(£-2w(g + -)) 



2 + 32ugq 
+ — (1 - ft) ^2gwfi - I w - fA ^. 



(g+-)ft-g 



.M i 

E 



»=i 



Z- Zi 



(4.17) 



The left hand side is a constant and the right hand side is a meromorphic function with 
simples poles at z = Zi and singularity at z = oo. The residues of the right hand side at 
the simple poles z — z% are 



Res, 



M 4 

p-^l^jj^i \ Z i ~ z p)\ z i ~ z l)\ z i ~ z j) 



-64g 2 



^-(n-l)z 2 + (q+l)z t 



M 

E 



^(Zi-Z^iZi-Zj) 



+ \ 4w 2 (ft 2 - 3ft + l)zf + IQujg 

M 2 



3(q+-)n-3q-l 



1 1 JVI 



jyi z * % 



ur 



+2— ft(l - ft),? 2 + 



8w 2 g(l - ft) + 8w 2 (g + -)(1 - Qf 



+4u{E - 2u{q + -)) 



^j + 32ugq 



(q + -)Q-q 



(4.18) 



Using the identities 

M M 



J2E- 



0. 



=1 jjti Z i Z j 



M M 

1=1 j^J 



-M(M-l), 



1 



i=l l^j^i \ Zi ~ Z l)\ Z i ~ Z j) 
M M 

E E 



Zi 



1 



( z i — z p){ z i — z l)( z i — z j) 
Zi 



i=\ p^l^j^i 
M M 

E E ■ 

i=l p^ljtjjti \ z i ~ z p)\ z i ~ z l)\ z i ~ z j) 



M M 

i=\ l^j^i \ Zi ~ z l)\ z i ~~ z j) 



0, 
0, 



0, 



(4.19) 



we can show that 



1, 



E-2u(q--)) -A 2 -Aoj 2 q 2 (l-n) 2 



M 



M 



E ^^ + 4w 2 (ft 2 - 3ft + l)M(M - 1) + 2— ft(l - ft) E<^ 



j=i 



2;- Zi 



!J 



+ 



8w 2 g(l - ft) + 8w 2 (g + -)(1 - ft) 2 + Auj(E - 2w(g + -)) 



i=l 



a; 



+ — (1-ft) 
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(2M + 2g)wft --u-E 



z. 



(4.20) 



The right hand side of (4.20) is a constant if and only if the coefficient of z as well as all 
the residues at the simple poles are vanishing. We thus obtain the energy eigenvalues, 



E 



~-w + 



2M - 2 ( q - I) + I 



wft, 



(4.21) 



and the constraint for the model parameters A, a; and g, 



A 2 + 4w 2 (l-ft) 



U! 



M 



M(M + 2q-l) + — ftE^ 

" i=l 



0. 



(4.22) 



Here the roots Zi satisfy the following system of algebraic equations, 

M 4 
2 2 V^ 

p^l^j^i \ Z i ~ Z p)\ Z i ~ z l)\ z i ~ z j) 



+ 9 



w(ft - l)z 2 + 4#( g + -)^ 



M 

E 



(Zi - Zi)(Zi - Zj) 



UJ 






ug 



1 



3(g + -)Q-3g-l 



Zj + 4# 2 g(g 



M 

E 



yWfi + (g+-)fi(fi-2) 



+ 2uigq 



(q+-)Q-q 



0. 



1,2, 



M 



(4.23) 



Here we have also used the relation (4.21) in obtaining (4.22) and (4.23). The correspond- 
ing wavefunction component if)+(z) °f the system is given by 

M 



M*) 



e 4 » 



(1-Cl)z 



i=l 



(4.24) 



and the component ip-{z) 



e A a 



as(i-ft)* 



(f)-(z) with (f)-(z) computed from (4.12) for A^O, 



4>-{z) = —j^C2-p4>+(. z )- Thanks to the exact solvability of £2- p , <f>-(z) automatically 
belongs to the same invariant subspace as (f)+(z). 

Some remarks are in order. The results obtained in [25] on the solution of the 2-photon 
Rabi model are incorrect. The reason is simple. As the author himself noted in the paper, 
the solution for the 2nd component of the two-component wavefunction of the 2-photon 
Rabi model does not belong to the same invariant subspace as the 1st component. There- 
fore his two-component wavefunction does not satisfy the matrix Schrodinger differential 
equation of the model and thus is not a solution to the coupled equations. 

As an example to the above general expressions, let us consider the Ai — 1 case. The 
energy is 



E 



1 



-co 



(2q + 2)cutt. 



(4.23) becomes 



to + q+(q + -)£l(n-2) 



Zi + 16 g q 



(q + -)n-q 



The solutions to this equation are 



«1 



4gq 



4gq{l - Q) - 2gfl 



wfi w(l-fi) 

Substituting into (4.22) gives the constraints A = and 

8 (g+£)^ 2 , 



(4.25) 

= 0. 
(4.26) 

(4.27) 
(4.28) 



A 2 + 8qco 2 

respectively. The constraint A = corresponds to the case of degenerate atomic levels. 
The constraint (4.28) agrees with those obtained in [19] by the Bogoliubov transformation 
method. The corresponding wavefunction ip+{z) is given by 

4^(1 -ft) -2gVt" 



M*) 



-£-(l-Cl)z I 



w(i-n) 



(4.29) 
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5 Quasi-exact solvability of the two-mode quantum 
Rabi model 

The Hamiltonian of the two-mode quantum Rabi model reads 

H = cu(a\ai + 0^02) + A<r 2 + ga x (a\al + 0102), (5.1) 

where we assume that the boson modes are degenerate with the same frequency u. In- 
troduce the operators K±, K , 

K + = alal, K_ = ai<22, K = -{a\ai + a\a2 + 1). (5.2) 

Then the Hamiltonian (5.1) becomes 

H = 2u (k - i) + Aa z + ga x (K_, + K_). (5.3) 

The operators K±,K form the usual su(l, 1) algebra. The quadratic Casimir of the 
algebra, C = K + K_ — K (K — 1), takes the particular values C = k(1 — k) in the rep- 
resentation (5.2), where k = 1/2, 1,3/2, • • •. This is the well-known infinite-dimensional 
unitary irreducible representation of su(l, 1) known as the positive discrete series T> + (k). 
Thus the Fock-Hilbert space decomposes into the direct sum of infinite number of sub- 
spaces % K labeled by k = 1/2, 1, 3/2, • • •. 

The basis state in the subspace "H K , denoted as \k, n), n = 0, 1, 2, • • •, has the form 



,t\n+2K-l/ /7 t\n 

k»)= ^ { Z , 10), (5-4) 

(n + 2n- l)\n\ 



and the action of {K±, K Q } in this representation is given by 

K \k, n) = {n + k)\k, n), 

K + \k,u) = y/(n + 2K)(n + l)\K,n + l), 

K_\K,n) = y(n + 2k — 1)ti\k, n — 1). (5.5) 

Using the Fock-Bargmann correspondence (4.6), one can show that the infinite set of 

monomials 

z n 
VkM= i n = 0,1,2, •••, (5.6) 

(n + 2«-l)!n! 



forms the basis in the Bargmann subspace associated with the representation (5.5). Thus 
the operators K±, K (5.2) have the single- variable differential realization, 

A A^ A 

K = z— + k, K + = z, K_ = z— + 2k—, k = 1/2,1,3/2,-... (5.7) 

az dz A az 
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By means of the differential representation (5.7), we can express the Hamiltonian (5.3) 
as the 2nd-order matrix differential operator 



/f=2w (4 +K 4) +Aff = +9 ^( z+ ^ +2 4)' 



(5.i 



Working in a representation defined by a x diagonal and in terms of the two-component 

(ib (z) \ 
, we see that the time-independent Schrodinger equation 
i)_(z) j 

yields the two coupled differential equations, 



d 2 d 

d 2 d 

9 Z -T2^-( Z ) + 2(-wz + gK)—ip-(z) + 



gz + 2u I K 



E 



iP+(z)+Aip.(z) = 0, 



gz — 2u ( k 



E 



ip_(z) -Aip + (z) = 0. 
(5.9) 



If A = these two equations decouple and reduce to the differential equations of two 
uncoupled two-mode squeezed harmonic oscillators which can be exactly solved separately 
[24]. For this reason in the following we will concentrate on the A^O case. 
With the substitution 



ip±(z) = e s (1 A)z (f±{z) 



A 



'-& 



(5.10) 



where - < 1, it follows, 

f d 2 . t . d . 1 

< gz—— + 2\uiAz + ok - — \- 2ku>A — u — E } ip + = — Ao2_, 

[ dz 2 dz J 

r $ d 4u 2 1 

gz— + 2[w(A - 2)z + gn} — H (1 - A)z + 2kw(A - 2) + w + ,5 }> </?_ ■--= Ay? + . 



dz 2 



Then the differential operator 



dz g 



d 2 



d 



Li-m = gz—rr, + 2[o;Az + gn\- — h 2kloA — co — E 
dz 2 dz 



(5.H) 



(5.12) 



in the 1st equation is exactly solvable. Eliminating y?-(z) from the system, we obtain the 
uncoupled differential equation for <p+(z) 

oz— + 2[w(A - 2)2 + ok]— + (1 - A)z + 2ku(A - 2) + w + E \ 



dz 2 



dz g 



id 2 d 1 

x I gz — + 2[coAz + gn]— + 2kwA - u> - E \ ip+(z) = -A 2 ip + (z). (5.13) 



12 



This is a 4th-order differential equation of Fuchs' type. Explicitly, 



o 2 & <P+ , a 2 



9 z 



dz 4 



+v 









+ J4w 2 (A 2 - 3A + l)z 2 + Aug 



u 



+ 8— A(l-A)z 2 + 



3(k+-)A-3«- 1 
1, 



z + 4g k(k 



U] d?V+ 

r] dz 2 



Wk{\ - A) + 8w 2 (/t +-)(!- A) 2 



-4oj(E — 2uk)] z + 8ugn 



[k + -)A-k 



\ d<f+ 



. w 



+ ^ 4— (1 - A) (2kuA -u-E)z 
{ 9 



+4wV(l-A) 2 - £-2w(«--) +A 2 \<p 



0. 



(5.14) 



As we show below, this equation is quasi-exactly solvable provided that the system pa- 
rameters A,u and g satisfy certain constraints, and exact solutions are polynomials in 
z. To this end, we seek polynomial solutions of the the form to the differential equation 

(5.14), 



M 



<P+(z) 



H(z-Zi 

1=1 



M = 1,2, 



(5.15) 



where Ai is the degree of the polynomial and Zi are roots of the polynomial to be deter- 
mined. Following the procedure similar to that in the last section, we obtain the energies 



1 



2M + 2 ( k - - ) + I 



E = -u + 
and the constraint for the model parameters A, u and g, 

A 2 + 4w 2 (l-A) 



uA, 



2u M 
M(M + 2k - I) + —AY Zi 
9 U. . 



0. 



(5.16) 



(5.17) 



Here the roots z% satisfy the following system of algebraic equations, 
M 4 

2 2 V"^ 

p^l^j^i \ Z i ~~ z p)\ z i ~~ z l)\ z i ~ z j) 



+ 4g 



1. 



oo{A-l)zf + g(K+-)z 



M 

£ 



(Z { - Zi)(Zi- Zj) 



+ i 4w 2 (A 2 - 3A + l)zf + Aug 



3(k + -)A-3«- 1 



+ V«(«+o)rE 



M 



2 J -^ Zi Zj 



U! 



+ 8— A(l - A)zf 
9 



■Sid 



MA + (k+-)A(A-2)+k 



Zi 



+ 8uigK 



[k + -)A-k 



i = 1,2, •••At. 



(5.11 
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Here we have also used the relation (5.16) in obtaining (5.17) and (5.18). The correspond- 
ing wavefunction component ip+{z) is given by 



ip+(z) = e » 



(1-A)z 



M 



H(z-Zi 
i=l 



(5.19) 



and the other component is ip-(z) = e~~° ~ ' z <p^(z) with (f-(z) determined from the first 
equation of (5.11) for A^O, V-0 2 ) = ~~ \£-i-m ^+( z ). Since L<i-m is an exactly solvable 
operator, tf-(z) automatically belongs to the same invariant subspace as tp + (z). 

As an example of the above general expressions, we consider the M. — 1 case. The 
energy is 



E = -w+(2k + 2)uA. 



(5.20) 



(5.18) becomes 












w 2 A(l - A)z\ + ug 


A + k + (k+-)A(A-2) 


Zi + g 2 n 


"(« + ^)A- 


- K 


= 0. 






(5.21 


It has two solutions 


Kg qk{\ - A 


A) ' 






(5.22 



Substituting into (5.17) gives the constraints A = and 



1 



A 2 + 8klu 2 = 8 ( k + - ) oj 2 A 2 = 0, 



(5.23) 



respectively. The constraint A = corresponds to the case of degenerate atomic levels. 
The constraint (5.23) is the non-trivial one. The corresponding wavefunction ip+(z) is 
given by 

Mz) = e-^- A) '(z- 9K{1 - A) - 9A ). (5.24) 



w(l-A) 
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